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Abstract . In the recent literature on dark energy (DE) model building we have learnt 
that cosmologies with variable cosmological parameters can mimic more traditional DE 
pictures exclusively based on scalar fields (e.g. quintessence and phantom). In a pre- 
vious work we have illustrated this situation within the context of a renormalization 
group running cosmological term, A. Here we analyze the possibility that both the 
cosmological term and the gravitational coupling, G, are running parameters within a 
more general framework (a variant of the so-called "AXCDM models") in which the 
DE fluid can be a mixture of a running A and another dynamical entity X (the "cos- 
mon") which may behave quintessence-like or phantom-like. We compute the effective 
EOS parameter, coe, of this composite fluid and show that the AXCDM can mimic 
to a large extent the standard ACDM model while retaining features hinting at its 
potential composite nature (such as the smooth crossing of the cosmological constant 
boundary LVg = —1). We further argue that the AXCDM models can cure the cosmo- 
logical coincidence problem. All in all we suggest that future experimental studies on 
precision cosmology should take seriously the possibility that the DE fluid can be a 
composite medium whose dynamical features are partially caused and renormalized by 
the quantum running of the cosmological parameters. 
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1 Introduction 



One of the most alluring aspects of modern cosmology is that it has become an accurately testable 
phenomenological branch of Physics. This has granted the field a fairly respectable status of 
empirical science, which certainly did not possess some two decades ago. Undoubtedly the most 
prominent accomplishment of cosmology to date has been to provide strong indirect support for the 
existence of both dark matter (DM) and dark energy (DE) from independent data sets derived from 
the observation of distant supernovae [1] , the anisotropies of the CMB [2] , the lensing effects on the 
propagation of light through weak gravitational fields [3] , and the inventory of cosmic matter from 
the large scale structures (LSS) of the Universe [Ij. But in spite of these outstanding achievements, 
modern cosmology still fails to fulfill the most important of its tasks, to wit: unraveling the 
ultimate physical nature of the components that build up the mysterious dark side of the Universe. 
Ignoring for the moment the puzzling DM sector (which is recently receiving spectacular empirical 
support [3) and whose final identity might hopefully be elucidated by directly detecting some day 
an exotic form (e.g. supersymmetric) of elementary particles, the DE component is nonetheless 
the most intriguing, unnerving and distressing of all. It turns out to be constant or slowly varying 
with time, it does not cluster, it fills smoothly all corners of our patch of accessible Universe 
and - to pitch even higher our level of perplexity - it constitutes by far the dominant form of 
energy density at the cosmological scale. What is it? Is it vacuum energy in the form of a 
cosmological constant? Is it the value of a slowly evolving homogeneous and isotropic scalar field? 
Perhaps a remnant of higher order gravitational terms in the effective action? Or some hint of 
a modified gravitational theory?... We don't know. But we do know at least (as of the early 
work of Zeldovich in 1967 [6]) that the typical size of the vacuum energy in QFT is exceedingly 
large to be reconciled with the characteristic energy densities at the cosmic scale, such as the 
cosmological matter-radiation density, the critical density, and of course also the cosmological 
energy density [p\ = A/SvrG) associated to A. All of them are of order of (10~^eF)^, far too 
small compared to the average energy densities in particle physics, say in the standard model of 
strong and electroweak interactions. Unfortunately, moving from QFT to string theory does not 
seem to help much, for after the process of compactification from 11 dimensions down to 3 we are 
left with a vastly complex "landscape" consisting of some 10^'''^'^ metastable (non-supersymmetric) 
vacua where to entertain our choice of the ground state [7] . All in all it looks very hard to identify 
p\ with a vacuum energy density. If we insist on it, we stumble once again upon the excruciating 
cosmological constant (CC) problem whose only known "technical" solution at present (though 
certainly not a very natural one!) is to postulate the existence of an ad hoc series of extremely 
fine-tuned cancellations among the various contributions [HI [9] §. 

In view of this situation it is perhaps wiser to avoid subscribing for the moment to a too strong 
opinion on the ultimate nature of the A term in Einstein's equation. Instead, we may treat it 
phenomenologically as a parameter in QFT in curved space-time, therefore acquiring the status 
of an effective charge similar to any running quantity in, say, QED or QCD. If one takes this 
practical point of view [10], we may get some guidance from the powerful renormalization group 

■^See e.g. [10] for a summarized presentation. 
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(RG) methods and use them to extract testable information on the possible evolution of A [llj 
from its presently measured value - see e.g. Ref. [12] for a devoted attempt in this direction. 
The RG evolution of the CC is not primarily a time evolution, it stems rather from it being a 
function A = A(/i) of the typical energy scale of the Universe at any given stage of the cosmic 
evolution. Following '[iO\ we identify this energy scale with the expansion rate or Hubble function, 
= H, which is of course time dependent. Therefore A = A{H) becomes also time dependent. 
The resulting model is a FLRW type cosmological model with running cosmological parameters. 
This point of view has been amply exploited in different ways in references |10j-|15] and constitutes 
an alternate candidate to quintessence [16j. Both pictures (quintessence and RG cosmology) have 
in common the idea of a dynamical DE. This can be useful in that a time-evolving DE may help 
to understand another aspect of the CC problem which is also rather intriguing, the so-called 
"coincidence problem" ^ or the problem of why the presently measured value of the CC (or DE) 
is so close to the matter density, i.e. both being at present of order of (10~'^ eV)'^. This problem has 
been addressed from different perspectives within modified quintessence models |17] . Remarkably, 
RG models can also tackle efficiently the coincidence problem [18], a point that will be reexamined 
here too. 

One may describe RG cosmological models from the viewpoint of an effective scalar field 
parametrization of the DE and compute e.g. the effective equation of state (EOS) for the running 
A. This was done in detail in [19j for a particular case with fixed gravitational coupling G and with 
non-conserved CC and matter densities. In a subsequent work a general algorithm was established 
to compute the effective EOS for general cosmological models with variable cosmological param- 
eters |20) . A further generalization (the so-called AXCDM models ^Sj) is to admit the possibility 
that the DE is a composite medium made out of a running A and a dynamical component X that 
we called the "cosmon" in [isjfl. In the present work we elaborate further along this line but we 
allow the gravitational coupling G being variable. In this way we obtain a picture in which matter 
and cosmon densities are separately conserved, whereas A evolves together with the gravitational 
coupling G according to the Bianchi identity - thus preserving general covariance. Apart from 
A = A(H) we explicitly obtain the function G = G{H). The outcome is a variant of the AXCDM 
model in which the change of A is compensated for by the evolution of G. 



2 Cosmon models with running A and G 

Consider Einstein's equations with a CC term. A, and move it onto the r.h.s of the field equations 
to form the combined quantity 8 vr G T^u . Here T^i, is the effective energy- momentum tensor T^i, = 
T^u + g^u Pk, where T^y is the ordinary contribution from matter-radiation, and /3a = A/(87rG). 
Next let us impose the Bianchi identity, (^T^i^ = 0. Upon evaluating it in the FLRW metric 
one finds 

^ [G (pm + pa)] + Gi?Qm/Om = 0, = 3(1 + w^) , (2.1) 

■^The name "cosmon" was originally introduced in ,21, to represent a dynamically adjusting DE field. Here 
we borrow the name to mean any (eflFective or fundamental) dynamical component of the DE (not necessarily 
quintessence-like) other than the cosmological term. 
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where pm and /9a are the matter-radiation and CC densities, and uJm = 0,1/3 (a^ = 3,4) for 
matter and radiation respectively. Notice that we contemplate p\ = dpp^/dt ^ and moreover we 
do not drop G from (j2.ip because we also admit the possibility that G 7^ 0. Let us next assume 
that the DE fluid is a composite medium consisting of the CC term and of another dynamical 
quantity X ("the cosmon"), with energy density px, characterized by an EOS parameter ivx- 
Then it is easy to see that Eq. (12. Ih generalizes into 

^[G {pm + Pd)]+GH {amPm + axPx) = ax = 3(l+a;x), (2.2) 



where pn = PA + Px '^^ the total DE density. Following ^8] we call the resulting composite DE 
model the AXCDM model. There are, however, different possible implementations of the AXCDM 
model: if we assume that pD and pm are separately conserved, then G must be constant and we 
obtain precisely the sort of AXCDM framework studied in detail in Ref. [18] , which will be referred 
hereafter as type I cosmon models. Alternatively, if we take as conserved quantities px and pm, 
then we must have pi + ai piH = for both pi = pm , px , and in this case ()2.2p implies that G 
must vary in time in combination with p/^ according to the equation 

G{pA + PX + Pm) + GpA = 0. (2.3) 

The class of AXCDM models based on these equations will be referred to as type II cosmon 
models and its study is the main aim of this work. We remark that the type II models constitute a 
generalization of the G-running cosmological model presented in Ref. [13j where the new ingredient 
here is the introduction of the self-conserved cosmon entity X which, as we will see, can play an 
important distinctive role. We should stress from the beginning that the cosmon need not be a 
fundamental field. In particular, X is not necessarily a quintessence scalar field |16]. not even a 
"standard" phantom field [22]. In fact, the barotropic index ujx of X can be above or below — 1, 
and moreover the energy density px can be positive or negative. Thus e.g. if uJx < — 1 and px < 
we obtain a kind of unusual fluid which has been called "phantom matter" (see Fig. 1 of [18j). 
Phantom matter does preserve - in contrast to usual phantom energy ~ the strong energy condition, 
as ordinary matter does. Clearly the cosmon plays the role of a very general entity ranging from 
the overall representation of a mixture of dynamical fields of various sorts, to the effective behavior 
of higher order curvature terms in the effective action. In spite of the phenomenological nature 
of this approach at the present stage, we cannot exclude the existence of a consistent Lagrangian 
formulation for some of these models. In the remainder of this paper we focus on the interesting 
phenomenological implications of the type II models. Still, it is important to say that the type II 
models are motivated by the models with a well defined Lagrangian formulation such as QFT in 
curved space-time [TU], to be discussed in the following paragraphs. 

To solve Eq. (j2.3p for G we need some input on the evolution of p\, otherwise we are led to the 
trivial solution G = const. Following \10\ [T2] we shall adopt the RG inspired evolution equation 



. M'pp\ (2.4) 
amp Air 

where is a free parameter: it essentially provides the ratio squared of the heavy masses contribut- 
ing to the /3-function of A versus the Planck mass, Mp. We naturally expect <C 1. As mentioned 
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in the introduction, we will choose fi = H as the typical RG scale in cosmology [lUj . Assuming a 
FLRW type of Universe with flat geometry we finally meet the following set of equations for the 
type II cosmon models: 

2 8tt G . , 

H = {Pm + PA + PX) , 

PA = Co + Ci , 

{Pni + PA + Px)dG + GdpA = 0. (2.5) 

The first one is Friedmann's equation for a spatially flat Universe, the second equation is the 
integrated form of (j2.4p with the boundary condition p\ = (the measured value of the CC 
density) for H = Hq (Hubble's parameter at present), hence 

Co = pl-^Mj,Hl C, = ^^Ml- (2.6) 

and the third equation is an equivalent differential form of ()2.3p . The basic set (j2.5p can be 
analytically solved to determine G as a function of the scale p = H, with the following result: 

G{H) = — , (2.7) 
1 + 1/ In ^ 

where Go = 1 /Mp . This logarithmic running law for G formally coincides with the one obtained in 
the framework of [13]. Of course the time/redshift dependence of G will be different here because 
of the cosmon contribution. The form of (j2.7p suggests that acts also as the /3-function for the 
RG running of G. To obtain G = G{z) and p\ = pa{z) as functions of the cosmological redshift z 
requires more work on solving differential equations. We limit ourselves to quote the final results. 
The gravitational coupling can only be expressed as an implicit function of z: 

1 . / 1 



9{z) 



1 + uln (^-ly =uln [nl, (1 + z)"™ zT^ ^Vl\-v\, (2. 



where we have defined (7(2;) = G{z)/Gq and $7^ = p^/p^ for the various components, with the 
critical density at present. Equation (j2.8p defines implicitly the function g = g{z), and once this 
is known the CC density as a function of the redshift is obtained from 

/ ^ Pa + J^ {pmjz) + pxjz)) g{z) -vpl 
Pa{z) = ^-^^ . (2.9) 

With these equations the type II class of cosmon models is solved. However, a highly convenient 
next step to do is to characterize this model with an effective equation of state (EOS) for its DE 
density. This is most useful in order to better compare the effective behavior of this model with 
alternative models of the DE (e.g. quintessence type models). 



3 Effective equation of state for the DE in type II cosmon models 

The effective EOS for type I cosmon models was studied in great detail in [18]. Let us thus 
concentrate here on type II models only. The general procedure to obtain the effective EOS 



5 



parameter cjg for models with variable cosmological parameters was thoroughly explained in |19| 
and [20] and we refer the reader to these and other useful references on this subject [23] . The 
formula for uje as a function of the cosmological redshift is 

UJe{z) = -1 + —— —— — . 3.1) 

3 Pd{z) dz 

The density po (not to be confused with /9_d) is the total DE density in the effective DE picture |20] . 
The relation between p^ and po can be obtained from matching (i.e. equating) the expansion 
rates H in the two pictures (in this case the AXCDM model and the effective DE picture): 

G{pm + Pd) = Gq {pm + Pd) ■ (3.2) 

By definition, the effective DE picture [20j is characterized by self-conserved DE and matter densi- 
ties {pD, Pm) and by a constant gravitational coupling: G = Gq. Generalizing the procedure of [20] 
we can show that pd{z) satisfies the following differential equation (stemming from the Bianchi 
identity) 

dpD Pd{z) -i{z) 

—r- = OLm — , (3.3) 

dz I + z 

where 

e(z) = giz) r pa(z) + "'"""^ p'x (1 + . (3.4) 

\ am J 

The solution of (j3.3p satisfying /Oz)(0) = p\, reads 

dz' i{z') 



PD{z) = {\-rZr 



Pd - Or 



(3.5) 



/o (1 + z')^°™+^^ 

Again following the methods of |20j it is convenient to rewrite Eq. (13.21) as follows: 

Mz) (1 + zr- + fA{z) + fx{z) (1 + zr^ = 0^ (1 + zr- + noiz) . (3.6) 

Here CId(z) = pd{z)/p^ with pd{z) given by (|3.5p . The fi functions above, whatever it be their 
form, must satisfy fm{0) = /a(0) = /x(0) = 1 in order to preserve the cosmic sum rule of the 
AXCDM models, which reads (in the flat case) 

n'i + n% = Q'^,, + nl + n^x = l. (3.7) 

At the same time the parameters (l^ and (l^ in the effective DE picture fulfill their own cosmic 
sum rule Q'^^ + n% = I. In general the quantities AJ7^ = 0^ - An% = n% - ^1% will be 
non- vanishing because they correspond to different parametrizations of the same data \19\ [20] . Of 
course it only makes sense to distinguish between Q,^ and (l^ when we fit densities of matter and 
of DE e.g. from distant supernovae data. However, for the radiation component (which is very 
well determined by CMB measurements) we set these differences to zero (see Section [3]). 

We may now compute dpo{z) / dz from (j3.6p and insert the result in (j3.ip . In the process we 
use the Bianchi identity derived in the previous section and the conservation laws in the AXCDM 
models. For the type II cosmon model it is easy to see that these conservation laws actually entail 
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fmiz) = fxi^) = 9{z)- And moreover the Bianchi identiy insures that the fohowing differential 
expression vanishes identically: 

nlfji + zT- + nl^ + n'^fji + = o . (3.8) 

Using these relations we can obtain, after a straightforward calculation, the final result. Let us 
present it in compact form as follows: 

u;,{z) = -1 + -5^ (3.9) 

where 

6iz) = (^g{z) - 17^) (1 + z)^- + ax g{z) ^\ (1 + z)^^ . (3.10) 

As expected, for the special case characterized by = (no running A), = (no cosmon) and 
AJl^ = (no difference between the parameters of the two pictures) we obtain 5{z) = 0, for all z, 
and we retrieve the EOS parameter of the CC: tOe = —1- Recall that the function g{z) = G{z)/Gq 
has been obtained implicitly in Eq. (|2.8p and that it takes the simple explicit form ()2.7[) only 
when written in terms of H. The latter (log) form suggests that G{z) evolves very little with z and 
therefore we expect that the departure of the effective EOS parameter ()3.9p from the CC boundary 
will be mainly controlled by the parametrization difference AfiJ^ and by the cosmon contribution - 
the last term on the r.h.s. of (j3.10p . We shall illustrate the behavior of (j3.9p with some non-trivial 
numerical examples in Section[5l Let us now study the bounds on i/, the parameter that regulates 
the evolution of both p\{z) and G{z). 



4 Nucleosynthesis bound on u. Potential implications for cosmo- 
logy and astrophysics 

As mentioned above, the type II cosmon models furnish a logarithmic law G = G{H), Eq. (j2.7p . 
which is formally identical to the one obtained in Ref. [13] in which the cosmon was absent. In the 
last reference a bound was obtained on the parameter v by considering the experimental limits 
placed by nucleosynthesis on the variation of G, with the result |z^| < 10^^. Here we are going 
to obtain a more stringent bound on i/ from the experimental restriction on the ratio of total 
DE to matter-radiation densities in the effective DE picture f = pD/Pm = ^D/^m- This ratio 
evolves with cosmic time or redshift. Particularly, in the standard ACDM model this ratio is given 
by p\/ Pm and increases without end as time passes by because pm — 0. At present its value 
is ro = p\/ p%i — 7/3, i.e. it is now "coincidentally" of order 1. This is in essence the cosmic 
coincidence problem mentioned in the introduction: why is this ratio of 0{1) right now? Let us 
denote the value of r at the nucleosynthesis epoch by = {po / Pm)N ■ As in previous works 
we will use the upper bound \fN\ ^ 10% \1'6\ I24| . For higher values of the expansion rate 
at nucleosynthesis would be too large and the amount of primordial helium synthesized would 
overshoot the experimental limits. To convert the upper bound on [fjv| into an upper bound on 
we substitute ()2.7p in (j2.8p and solve the resulting equation for v. The final (exact) analytical 
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rGsult IS 

Let us use the formula for the expansion rate in the effective DE picture to ehminate H'^(z)/Hq 
from (j4.ip in terms of the ratio r defined above: 

H^{z) 



Hi 



n^{z)+nD{z) = nl,{l + zr-^{l + f). (4.2) 



At the nucleosynthesis epoch (zjv ~ 10 ) we have am = 4 because radiation dominates. Being 
z very large we can neglect il^ in the numerator of (j4.ip . Furthermore, given the fact that 
X is assumed to be a DE component {—1 — 5 < lux < —1/3, with 6 > small), we have 
ax = 3(1 + ujx) < 2 (hence ax < Om) and the cosmon contribution can also be neglected at 
z = zx- Finally, as mentioned in Section [3l we set AOj^ = for radiation. All in all at the 
nucleosynthesis epoch Eq. (|4.ip boils down to 

(4.3) 



-fx^\{l + Zx)^ -fx 



1 + nlil + zxY{l + fx) [In + zxY{l + fx)) - 1] (1 + fN) In [^^ z%] 

Here fJ^j is the total radiation density (from photons and neutrinos) normalized to the critical 
density at present; il^j is given by 1.68 times the standard value of the normalized photon density 
ilJJ ~ 4.6 X 10~^ (for h = 0.73) [2]. From these equations the desired nucleosynthesis bound on v 
ensues immediately: 

\v\<lQ-\ (4.4) 

which is an order of magnitude more restrictive than the previous bound obtained in [T3]. Since 
the bound (14. 3p essentially does not depend on the X (cosmon) component, it also applies to the 
models studied in [13]. The fact that u <^ 1 \s indeed a natural theoretical expectation from 
the interpretation of this parameter within effective field theory (cf. Ref. |12j for details). An 
obvious implication of this bound on type II cosmon models is that the corresponding running of 
A is in principle rather hampered as compared to that in type I models [15] . In the latter, the 
nucleosynthesis bound affected a combined parameter involving both the A and cosmon dynamics: 
e = ^{l + uJx) ^ 0.1. However, for type II cosmon models nucleosynthesis does not place any 
bound on the cosmon barotropic index lox and in this sense we have more freedom to play with 
it. We point out that the ratio r in the AXCDM picture, i.e. r = PD/Pm, is related to f as 
r = g{z)~^ {Q,^/n^) (1 + f) — 1, where we have used (|3.2p . Let us next consider this expression 
at the nucleosynthesis time, where H{zx) = Hx- Setting Afi^j = 0, as previously discussed, and 
using (j2.7p we find 

rx = fN + v{l + fx) lnH%/Hl ^fx + ^{1 + rx) In [n^p. z%\ . (4.5) 

In view of (|4.4p the second term on the r.h.s. of the above expression remains below 10%. Therefore 
we find that at the nucleosynthesis time the bound on r lies within the same order of magnitude 
as the bound on f, i.e. |rAr| < 10%. Taking this into account and making the natural assumption 
that the value of AO^ is small (big differences between the two pictures are not expected), we see 
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that in practice we need not distinguish between f and r. This is exphcitly corroborated in the 
numerical analysis performed in the next section^. 

The tighter upper limit (|4.4p on type II models from nucleosynthesis as compared to type I does 
not necessarily go to the detriment of the potentially important role played by u when considering 
other implications. For example, for sufficiently small v > Q (specifically for v ~ 10-6), the type 
II cosmon models can offer an explanation for the flat rotation curves of galaxies exactly as the 
RG model of Ref. [13j[f|. This is because the type II models preserve the running logarithmic law 
(|2.7p wherein at the astrophysics level H^^ is replaced by the radial coordinate of the galactic 
system [13] . Remarkably, we see that even for very small (but positive) values oi v - far below 
the limits placed by cosmology considerations (j4.4p - we may have dramatic implications at the 
astrophysical scale. 

For 1/ < 0, however, these astrophysical implications are not possible. But in compensation 
we can get a handle on the cosmological coincidence problem mentioned in the introduction. The 
line of argumentation closely follows the approach of Ref. [18] for type I models. Namely, if the 
Universe's evolution has a stopping (or, more properly, a turning) point at some time in its future, 
then the ratio r will be bounded from above for the entire lifetime of the Universe, and therefore 
our present time cannot be considered special; that is to say, in a AXCDM model r can be of 
order 1 not because - as argued within the standard ACDM model - we just happen to live at 
a time near the transition from deceleration into acceleration, but simply because r can never 
increase beyond a fixed number. Moreover, in a large portion of the parameter space r remains 
below 0{1 — 10). We refer the reader to [18J for an expanded exposition of this approach to the 
coincidence problem, here we limit ourselves to apply the same arguments within the context of 
type II models. For these models the turning point associated to the possible resolution of the 
coincidence problem exists only for v <Q. This can be seen immediately from ()2.7p : H will vanish 
in the future at the very same point where G also vanishes. However, the stopping redshift value 
(call it zs) can be approached with the correct sign of the gravitational coupling (G > 0), if and 
only if 1/ < 0. For > we would have G < before reaching zs- One can also show that the 
existence of the turning point can be realized both for the cosmon barotropic index ojx above or 
below —1, i.e. for quintessence-like or phantom-like cosmon. To study the precise conditions for 
the existence of the turning point it proves convenient to use Eq. (j4.ip and set H{z) = in it. The 
root z = Zs satisfying the resulting equation will obviously lie in the matter dominated epoch, so 
we can set am = 3 in (j4.ip . Introducing the (continuous) function 

fiz) = nlii + zf + + .)3(i+--) + (1 - 0^ - 0°, - z.) (4.6) 

we find, equivalently, that zs is obtained as a root of it: f{zs) = 0. Next we apply Bolzano's 
theorem for continuous functions, and to this end we explore the change of sign of f{z) at present 

''We point out that the bound on the relative abundance of DE in the nucleosynthesis epoch [24, applies originally 
to f because in that analysis it is assumed that G is constant. Since, however, this implies a tight constraint on i/, 
Eq. (|4.4|l . it follows that the bound on r essentially applies to r as well, at least within order of magnitude. 

^ We omit details here. See section 7 of Ref. [13] for a detailed discussion of the possible astrophysical implications 
of a running G on the flat rotation curves of galaxies. Type II cosmon models lead to the same kind of picture. 
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{z = 0) as compared to the remote future {z = —1). We find 

n./,.) . ('-"'"-f-^ " (4.7) 

[ sign(SZ^) ■ oo if < — 1 
/(O) = 1-1/ (4.8) 

Being 1 — obviously positive (because z/ < 0) we can guarantee stopping point zs in the future 
under the following conditions: 

0^>l-0[^-i/ for ujx>-l, (4.9) 

n°x<0 for wx<-l- (4.10) 

Notice that — 1 < cox < — 1/3 corresponds to a quintessence-like cosmon. Given the approximate 
prior = 0.3 (e.g. from LSS [4J and supernovae data[l]) and using the bound (j4.4p and the 
sum rule (|3.7p . we learn that the corresponding stopping condition (|4.9p is essentially equivalent to 
require < 0. This should not be considered as a drawback; within our composite DE framework 
the generalized sum rule (j3.7p makes allowance for any sign of Q^, provided the total DE density 
$1^5 = + ^x positive and of the order of the experimentally measured value ~ 0.7. In 
other words, this essentially implies Qj^ > and it confirms that the condition (j4.9p is realized by 
a quintessence-like cosmon whose density is, therefore, decreasing with time. Being the running 
of A very mild it is not surprising that the stopping condition is almost equivalent to < 0. On 
the other hand, for < and —1 — 6 < lox < — 1 we must necessarily have il^^ < 0, Eq. (|4.10p . 
and in this case the cosmon behaves like "phantom matter" rather than as usual phantom energy 
(cf. Fig. 1 of [E]). Finally, let us consider the simplest case u = (i.e. cosmon models with fixed 
A and G). In this particular instance the total normalized DE density is the same as in the type 
I cosmon models for z/ = |18j : 

nDiz) = nl + n''j,{i + zf^^+'^^K (4.ii) 

For uJx > — 1 (resp. lox < — 1) the stopping point exists if il^ < (resp. fi^ < 0). This reflects 
the continuity of our analysis with respect to the variable v as this situation corresponds to the 
z/ — > limit of the z^ 7^ case studied above. Remarkably, for any value of v the existence of 
a turning point for ujx < — 1 does not correspond to X behaving as phantom DE but rather 
as "phantom matter" [LS^ . For this reason there is no "Big Rip" singularity [22] in this kind of 
models. To summarize, the existence of a turning point is guaranteed in many circumstances 
within the framework of AXCDM cosmologies, and it appears to be a very desirable ingredient 
to help quenching the severity of the cosmological coincidence problem and at the same time to 
eschew future cosmic singularities of the "Big Rip" type. In its absence (namely of a turning point) 
the Universe's expansion would be eternal and the ratio r = pn/Pm would tend to infinity as in 
the standard ACDM model. In the next section we present some numerical examples illustrating 
the various possibilities described above. 
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Figure 1: Some examples of the behavior of the effective equation of state function ()3.9p for the 
type II cosmon models: (a) ujx = —1.65, v = +0.001, = 0.67, AJ]^ = 0.01. In this case, the 
EOS presents an effective transition from quintessence to phantom regime in the recent past; (b) 
uJx = —0.85, V = —0.001, fi^ = 0.3, Afi^ = —0.01. Here the transition is just the other way 
around; (c) ujx = —0.95, v = 0.001, 17^ = 0.75, Ail^ = 0, the EOS mimics the behavior of a 
cosmological constant; (d) ujx = -1.15, v = -0.001, Afi^ = and ^\ = 0.4 (solid) or Vt\ = 0.8 
(dashed). In the last two examples the EOS parameter remains in the phantom or quintessence 
regime for all the redshifts attainable by present and scheduled supernovae experiments. 



5 Numerical analysis 

The results of the preceding sections are further illustrated in this section through particular 
numerical examples. We present plots of the evolution of the dark energy effective EOS parameter 
a;e, the Hubble parameter H and the ratio of the energy densities of dark energy and matter 
components r. These plots depict the interesting phenomena that the type II cosmon models 
exhibit in different parameter regimes. We first focus on the redshift dependence of the DE effective 
EOS parameter Ue, illustrated in the plots of Fig. [TJ Each of these plots demonstrates how type 
II models can result in the features of uje{z) which are consistent with the current observational 
data. The type II cosmon models can exhibit a very interesting phenomenon of the CC boundary 
crossing, both the crossing from the quintessence regime to the phantom regime with the expansion 
and the transition in the opposite direction. In Fig. [T^ we show an example of the realization 
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Figure 2: The effective equation of state parameter cUe = uJe{z) for the case ujx = —1.85, 
v = —0.001, fi^ = 0.75 and different values of AO^. We see that this last parameter crucially 
determines the evolution. All the curves present stopping of the expansion and a ratio between 
dark energy and matter radiation bounded from above. 

of the crossing from the quintessence regime to the phantom regime for the following values of 
parameters: tox = -1-65, u = +0.001, = 0.67 and AOjii = 0.01. An example of the transition 
in the opposite direction for parameter values lox = —0.85, u = —0.001, = 0.3, Ail^ = —0.01 
is depicted in Fig. [TJj. The plot in Fig. ^ shows that the AXCDM models may mimic the behavior 
of the cosmological constant. The values of the parameters in this case are ujx = —0.95, ly = 0.001, 
il^ = 0.75 and = 0. Finally, in the plot given in Fig. [T]i we give two examples in which 

the type II models exhibit either quintessence-like or phantom-like features for redshifts amenable 
to the SNIa observations. An example of the quintessence-like behavior (the dashed line in Fig. 
[TJl) is obtained for lox = -1-15, u = —0.001, AQ^ = and Q\ = 0.8, whereas an example of the 
phantom-like behavior (the solid line in Fig. [T]i) is realized for tox = -1-15, u = —0.001, AfiJ^ = 
and Q\ = 0.4. These examples illustrate the potential of the AXCDM models in explaining various 
types of dynamics for u)e{z). 

The dependence of the form of the function iOe{z) on the value of the parameter AJ7^ is given in 
Fig. [2] for the following values of the remaining parameters: wx = -1-85, u = —0.001, 0,^ = 0.75. 
From the figure it is clear that iOe{z) function changes considerably when AQ^ varies at a percent 
level, i.e. it is rather sensitive to the choice of AJ7^. All curves in this figure are characterized by 
the stopping of the expansion of the universe and with r{z) = pd{z) / Pm{z) bounded from above. 

The discussion presented so far indicates that the AXCDM models can easily provide the 
effective dark energy EOS in accordance with the observational data. This feature is certainly 
expected since dark energy contains multiple components and a correspondingly larger number of 
parameters. This possible drawback is, however, more than compensated, by the characteristics of 
the type II models which offer a robust solution to the coincidence problem, as already discussed 
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Figure 3: The Hubble function in units of its present value and depending on (a) redshift or (b) 
cosmic time (measured in Hubble time units, H^^), for the values of the parameters of Fig. [2] and 
the specific choice AilJ^ = —0.005. The stopping and subsequent reversal of the expansion takes 
place after a time longer than the present age of the Universe, to = 13.7 Gyr ~ 0.99Hq^ . 

in detail in [18| for the type I models. The class of type H models studied in the present paper 
also offers the solution of the coincidence problem. Namely, it is possible to select the model 
parameters so that the expansion of the universe stops at some future moment. An example of 
such a scenario is depicted in Fig. O where the evolution of the Hubble parameter is given in terms 
of (a) redshift and (b) cosmic time. The model parameters used are those from Fig. [2] together 
with A17^ = -0.005. The evolution of the universe of this type is further characterized by the 
specific form of the evolution of the ratio r. Namely, r{z) (r(t)) is bounded from above, i.e. it grows 
with the expansion up to some maximal value and then starts to decrease, changes sign at some 
future moment and attains the value -1 at the stopping point. It is especially interesting that for 
a considerable volume in the model parametric space {ujx, ^^a) ^'^^ ratio r/ro (where denotes 
the present value of the ratio r) remains of the order 1 once the universe starts accelerating. This 
fact implies that the class of type H models can provide a solution (or a significant alleviation) of 
the coincidence problem. An example of the described dynamics of r in terms of (a) redshift and 
(b) cosmic time is given in Fig. [5] for the same set of parameters as in Fig. [sH . A very similar 
conclusion about the evolution of the ratio r was reached within the class of type I models studied 
in [IB]. The fact that both type I and type H models provide the solution to the coincidence 
problem indicates that this solution is not model dependent and gives support to a general claim 
that AXCDM models provide a robust explanation of the coincidence problem. 

Before finishing this section the following observation is in order. As we have seen, the type H 
models with flat FLRW metric and a prior on il^) (typically il^j = 0.7, the current value of the 
DE density) contain three parameters (wx, ^a)' J^^* ^ tyP^ I models [18]. Of these three 
parameters v is tightly constrained by the nucleosynthesis considerations, as described in Section 
m This parameter space can be compared with that of other DE models which also try to solve the 
coincidence problem. For instance, in interactive quintessence models (IQE) [T7| one has a similar 

^For a sufficiently small value of Af2^ (which is indeed expected to be so) and v values allowed by (|4.4|l . the 
plots for r and f are practically indistinguishable. For this reason the ones for r are not included in Fig|4l 
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Figure 4: The ratio f between dark energy and matter-radiation in units of its present value and 
depending on redshift (a) or cosmic time (b), for the values of the parameters of Fig. [2] and 
the specific choice Af][^ = —0.005. We see that the ratio presents a maximum that can help to 
explain or significantly alleviate the cosmic coincidence problem. The plots for r result to be nearly 
indistinguishable from these due to the small value of AQ^, and thus are not shown. 

number of parameters, to wit: the EOS parameter of the quintessence field, u^, the coupling of 
matter to quintessence (usually in the form of a source function - let us call HQ- which one 
has to introduce totally ad hoc and depends on at least one parameter), and finally one or more 
parameters related to the (assumed) form of the potential. Moreover, if one wishes to trigger a 
transition between quintessence and phantom regimes one usually has to resort to more complex, 
e.g. hybrid, structures [T7] in which at least one additional scalar field is necessary, though carrying 
a "wrong sign" kinetic term (i.e. it must be a ghost field). Similar considerations apply to A;-essence 
models [25] . where an ad hoc non-linear function of the kinetic terms must be introduced along 
with additional parameters. In our case we can modulate that transition through the parameter 
Ail^, as we have seen above. However, even this parameter is not a new input of the model since 
it actually belongs to the corresponding effective EOS picture, not to the AXCDM model itself. 
To summarize, quintessence-like models have in general a similar number of parameters (if not 
more) as compared to the AXCDM, and may present additional features that one has to ponder 
carefully. For example, in IQE models the aforesaid coupling source Q (whose presence is essential 
to tackle the coincidence problem) necessarily entails the non-conservation of matter. This is in 
contradistinction to the AXCDM models, both of type I and II, where matter is strictly conserved. 

6 Conclusions 

We have studied the possibility that the dark energy (DE) of the Universe is a composite medium 
where certain dynamical components could be entangled with running cosmological parameters. 
We have modeled this composite structure by assuming that the DE is a mixture of a running 
cosmological term A and a dynamical entity X (the cosmon). The latter is not necessarily an 
elementary scalar field, it could rather be the effective behavior of a multicomponent field system 
or even the result of higher order terms in the effective action. In fact the generality of X is such 
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that its effective barotropic index ujx can be both quintessence-Uke (wx ^ ~1) ^-iid phantom-hke 
[iox ^ —1), and the energy density of X can be positive or negative {px ^ 0). The kind of 
composite DE model we have studied is a variant of the previously considered AXCDM model |18j . 
In all of these models the total DE density, pu, is the sum of the cosmon density, px, and the 
A density, p\. However, in contradistinction to [IB], here the cosmon energy density is conserved 
(as the matter density itself), and thus the covariance of Einstein's equations (expressed by the 
Bianchi identity) requires that the running of A must be accompanied by a running gravitational 
coupling G. We have computed the effective equation of state (EOS) of this composite DE system 
and found that under suitable conditions it can mimic to a large extent the behavior of the 
standard ACDM model, but at the same time the fine details reveal the possibility to observe mild 
transitions from quintessence-like into phantom-like behavior, and vice versa. If observed in the 
next generation of high precision cosmology experiments (such as DES, SNAP and PLANCK [26|). 
it would point at the possible composite structure of the DE. Furthermore, we have revisited 
the cosmological coincidence problem within this model, and found that in a large portion of the 
parameter space (similarly to [18]) there is the possibility that the Universe exhibits a turning 
point in its evolution, a fact that would automatically keep the ratio between the DE density to 
matter density (r = pn/ Pm) within bounds. Incidentally, it would also free the cosmic evolution 
from future singularities such as the Big Rip [22] . We have displayed concrete examples where the 
ratio r stays within a few times its present value tq for the entire history of the Universe. This 
feature seems to be independent of the particular implementation of the AXCDM model (whether 
of type I or of type II as defined in our analysis) and it suggests that this kind of composite models 
could provide a clue to solving, or at least highly mitigating, the cosmological coincidence problem. 

Finally, let us mention that the choice of scale in the RG motivated approach depends on the 
symmetries of the underlying space-time metric. As an example, in reference [13] it was explored 
the possibility to apply the RG approach at the astrophysical level, which leads the RG scale to 
develop a dependence on the geometric coordinates of the system. Furthermore, when we move to 
the cosmological context we expect that the A term in an inhomogeneous Universe should develop 
spatial fluctuations which could affect the growth of the matter structure and even the CMB 
anisotropics. The detailed estimate of the effect of the A inhomogeneities on these observables 
requires an analysis far beyond the scope of this paper. To get a hint of the kind of investigation 
required we refer the reader to the work of [27J, which concentrates on a simpler RG model for A 
in which G is constant and there is no cosmon. In this case one finds that the RG running of A 
may have a significant effect on the structure growth, i.e. that the large scale structure data may 
provide strong constraints on the RG parameters. Although the RG effects in AXCDM models 
of type II are not expected to be so important, owing to the existence of the cosmon component 
and the severe nucleosynthesis constraints, the LSS and CMB data are important tools for placing 
further constraints on these models. This kind of considerations will be addressed elsewhere. 

To summarize, the AXCDM models seem to offer a viable extension of the standard ACDM 
model in which some cosmological problems seem to be better under control. At the same time 
they offer a new approach to modeling the DE behavior without compromising the analysis to a 
particular field structure of the theory (e.g. a specific scalar field potential). The dynamics of 
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these models is entirely determined by the RG behavior of the cosmological parameters together 
with the covariant conservation laws involving the different components of the DE. In our opinion 
the AXCDM models should be considered as serious candidates for describing the potential DE 
features observed in the next generation of experiments. 
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